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Abstract 

AMS-02 recently published its lepton spectra measurement. The results show 
that the positron fraction no longer increases above ~200 GeV. The aim of this 
work is to investigate the possibility that the excess of positron fraction is due to 
pulsars. Nearby known pulsars from ATNF catalogue are considered as a possible 
primary positron source of the high energy positrons. We find that the pulsars with 
age T ~ (0.45 ~ 4.5) x 10 5 yr and distance d < 0.5 kpc can explain the behavior 
of positron fraction of AMS-02 in the range of high energy. We show that each of 
the four pulsars — Geminga, J1741-2054, Monogem and J0942-5552 — is able to be 
a single source satisfying all considered physical requirements. We also discuss the 
possibility that these high energy er are from multiple pulsars. The multiple pulsars 
contribution predicts a positron fraction with some structures at higher energies. 
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1 Introduction 

The positron fraction spectrum e + /(e + +e _ ) in the cosmic ray (CR) contains two 
components: secondary e ± produced by nuclei collision and primary e - . It is cur¬ 
rently believed that these two components, each of which will produce a diffused 
power low spectrum, predict a positron fraction which goes down with energy. How¬ 
ever, the latest results measured by the Alpha Magnetic Spectrometer (AMS-02) 
with high accuracy indicate that the positron fraction increases with energy above 
~8 GeV and does not increase with energy above ~200 GeV BIZI. This “increasing” 
behavior, which is also observed by the payload for antimatter matter exploration 
and light-nuclei astrophysics (PAMELA) [3H6] and the Fermi Large Area Telescope 
(Fermi-LAT) [101 [TTj, is not compatible with only diffused power low components. 
The “cutoff” behavior above 200 GeV , which can be well described by a common 
source term with an exponential cutoff parameter in the Eq.(l) of |lj, indicates that 
potential sources produce the exceed of electron and positron pairs. 

AMS-02 [H 2] is a state-of-the-art astroparticle detector installed on the In¬ 
ternational Space Station (ISS). It carries a Transition Radiation Detector (TRD) 
and a Electromagnetic Calorimeter (ECAL). These two sub-detectors provide in¬ 
dependent proton/lepton identification, which will achieve a much larger proton 
rejection power of AMS-02 compared with PAMELA which has only one Electro¬ 
magnetic Calorimeter for proton/lepton identification using the 3D shower shape 
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and Energy-Momentum match (E/P). Compared with Fermi-LAT, AMS-02 has a 
large magnet which can identify charge sign of the particle. Thus, the contamination 
of electrons (also called “charge confusion” in [T]) in the positron sample of AMS-02 
is much smaller that that of Fermi-LAT. For the reasons given above, there is much 
less proton or charge confusion contamination in AMS-02 measurement than that in 
PAMELA or Fermi-LAT. Here, we only interpret AMS-02 result due to the lack of 
knowledge of the contamination control in PAMELA and Fermi-LAT measurements. 

The AMS-02’s recent measurements of positron fraction [Tj, e + flux, e~ flux 
[2] and (e~+e + ) flux [I2]were published. The e~ flux contains three components: 
primary e~, secondary e~ and e~ from unknown sources. The e + flux contains only 
two components: e + from secondary production and primary e + from sources to be 
indentihed. To avoid the unnecessary uncertainty of primary e~, the e + flux seems 
to be an ideal spectrum to study extra sources. However, there is an acceptance 
uncertainty from the detector itself in the e + and e~ fluxes. This uncertainty in e + 
flux is strongly correlated with that in e~ flux [2], especially at high energies. The 
positron fraction can avoid this systematic uncertainty [Tj. For example, one can 
clearly see a drop at the last point (350 GeV ~ 500 GeV) in the positron fraction 
but cannot tell a drop at the last point (370 GeV ~ 500 GeV) in the e + flux due 
to its larger error bars. Therefore, positron fraction is used to study extra sources 
while e - flux is used to estimate the primary e~ which will affect the denominator 
of e + / (e + +e“). 

Recent studies have proposed some interpretations, such as dark matter annihi¬ 
lation or decay pT3Tf2Tj . supernova remnants (SNRs) [22H2T| . secondary production 
in the interstellar medium (ISM) [28] and pulsars [20, 21, 1291 143] . Cosmic ray flux 
data can also be together with other observations (like the dark matter relic density 
and the direct detection experimental results etc.) to give a combined constraint 
on dark matter models mm- Besides dark matter scenario, the others can pro¬ 
vide astrophysical explanations which do not require the existence of new particles. 
SNRs model, for instance in [2T] and [TBj, introduce some new mechanisms for the 
propagation model or special distributions of the primary sources. The “model- 
independent” approach from [28], sets an upper limit of the positron fraction by 
neglecting radiative losses of electrons and positron but does not indicate any ob¬ 
vious cutoff in the spectrum. Among them, the pulsar interpretation is one of the 
scenarios which predict a cutoff at a few hundred GeV in the positron fraction spec¬ 
trum and do not contridict other cosmic ray spectrums (eg. boron-to-carbon). The 
pioneering works on pulsar interpretaion of positron fraction have been performed 
by [33] ]35s 32] a few years ago. Combined analyses of the recent AMS-02 lepton 
data have been performed by [27] and [20], with a global fit on positron fraction [I], 
e + flux, e~ flux and (e“+e + ) flux. To avoid the over-estimation of the y 2 , however, 
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only two out of four spectrums should be used in the fit. As the reasons given by 
the previous paragraph, we only study positron fraction and e~ flux in this paper. 

A pulsar is widely regarded as a rotating neutron star with a strong magneto¬ 
sphere, which can accelerate electrons, which will induce an electromagnetic cascade 
through the emission of curvature radiation [47l [oil] . This leads to the production of 
high energy photons which eventually induces e + e“ pair production. This process 
produces the same amount of high energy e + and e ~, which can escape from the 
magnetosphere and propagate to the earth. There is a cutoff energy of the photons 
produced in a pulsar, which leads to a cutoff in the positron fraction. 

In this paper, DRAGON joTI (ool is used as a numerical tool to model the prop¬ 
agation environment, to tune the related parameters and to estimate the e ± back¬ 
ground. The authors of [53TI56] did a very complete work on three-dimensional 
cosmic-ray modeling. In the 3-D models, they pointed out the spiral arms have an 
effect on the propagation parameters. A 2-D model is used in this paper because 
we focus on the lepton spectra implication. Due to the energy loss of leptons, the 
effect of spiral arms on the high energy leptons is less important than that of the 
additional nearby sources contribution. ROOT is used to minimize y 2 to get the 
best fit results. We consider six nearby pulsars from ATNF catalogue [57J J58j as 
the possible extra single sources of the high energy positrons. We find only four, 
which are Geminga, J1741-2054, Monogem and J0942-5552, can survive from all 
considered physical requirements. We then discuss the possibility that these high 
energy are from multiple pulsars. The multiple pulsars contribution predicts a 
positron fraction with some structures at higher energies. 

The paper is organized as follows. Section [ 2 ] shows the way where e ± background 
is estimated. In Section [3j the properties of pulsars are described and the profile of 
e ± fluxes produced by a pulsar is derived. The interpretation of positron fraction 
with one single pulsar is discussed in Section [4] and the hypothesis about multiple 
pulsars interpretion is tested in Section [5j The conclusions are drawn in Section [6j 
I 11 Appendix [Aj the diffusion energy-loss equation for a burst-like source is solved 
with the spherically symmetric approximation. 

2 Propagation parameters and e ± background 

The Galacitc background of the lepton fluxes are considered as three main com¬ 
ponents, which are primary electrons from CR sources, secondary electrons and 
positrons from the interactions between the CR and the interstellar medium (ISM). 
The propagation of e ± in the Galaxy obeys the following Ginzburg and Syrovatskii’s 
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equation ofT), also in [52J GTT 



dp p 2 dp 


r\ O ' O 


(p - |v • Vc)fi 


( 1 ) 


j>i 


where p = \p\ is the particle momentum; f t (x, t,p) is the particle number density of 
a species i per unit momentum interval; v c is the convection velocity, /3 = v/c is the 
ratio of velocity to the speed of light; <r in is the total inelastic cross section onto the 
ISM gas, whose density is n gas ; (r 3l is the production cross section of the species i 
by the fragmentation of the species j (with j > i); and Qi(x, t,p) is the source term 
of species i, which can be thought to be steady Qi = Qi(x,p ) for background CR 
particles. 

The spatial diffusion coefficient D in the cylindrical coordinate system (r, z ) may 
be parameterized as [52, 01), EHJ 



( 2 ) 


where p = pc/(Ze) is dehned as the particle magnetic rigidity, z t is the scale height 
of the diffusion coefficient, L is the halo size, and 5 is the index of the power-law 
dependence of the diffusion coefficient on the rigidity. Dq is the normalization of 
the diffusion coefficient at the reference rigidity p 0 = 4 GV. Previous DRAGON 
papers [53R55] tested a few models with the exponential profile, i.e., the left formula 
in (J2|, which is more physical than the constant one, i.e., the right one. The effect 
of choosing different profiles on the electron and positron background is small if the 
parameters are properly set. In this paper, the constant profile is used in order to 
compute the pulsar profile in an analytical way, i.e., eq. 0 in Section [3j The 
function f{r) describes a possible radial dependence of D, and it can be taken to be 
unity for simplicity. 


The diffusion coefficient in momentum space D pp is related to the spatial diffusion 
coefficient D by [18, 61| [62, 65] 


(3) 


D D = _ _ 

pp 3<5(4-5 2 )(4-5)u; 


where va is the Alfven velocity, and S is the power-law index as given in (J2]) . w is 


the ratio of magnetohydrodynamic wave energy density to the magnetic held energy 
density, and it is usually taken to be 1. 

DRAGON [5TH53] is used to tune the propagation parameters according to the 
B/C ratio, which is sensitive to the parameters. The Markov Chain Monte Carlo 
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algorithm (MCMC, |67j ) is used to determine D 0 and 6. The priors are shown in 
Table, [l] The posterior distributions can be shown in a contour in the Dq and <5 
plane in Fig. [l] 


Table 1: The priors of Dq and 5 



start value 

minimum value 

maximum value 

.D 0 (xl0 28 cm 2 s x ) 

3 

2.5 

7.5 

5 

0.40 

0.20 

0.65 



D q (10 28 cmV 1 ) 


Figure 1: Contour in the Dq and 5 plane. The cross shows the best fit value while 
the three closed curves from inside to outside show the 68.3% C.L., 95.4% C.L. and 
99.7% C.L. respectively. 
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The parameters with their 68 % C.L. uncertainties from the fit are as follows: 

(D 0 1 „ _ = ( 6.20 ± 0.31) x 10 28 cmV 1 

u I po=4 GV v ' 

5 = 0.31 ± 0.03 

(4) 

L = 4 kpc 7 

va = 40 krn/s 

where the halo size L is taken from the MED model of (63] . and the va is fixed. 
These parameters are consistent with what the authors of Ref. [2D] has got in the 
reaccelaration propagation model. 

To avoid the uncertainty of solar modulation, AMS-02 proton flux [3] above 45 
GV is fitted to get the injection spectra using MCMC [67]. Three breaks, which 
are 6.7 GV, 11 GV and 316(±148) GV, are introduced in the injection spectrum 
of nuclei. The proton spectral indice below and above the breaks are 2.25, 2.35, 
2.501 (±0.010) and 2.501-0.084(±0.050), respectively. The high energy spectral in¬ 
dices of helium, carbon and oxygen are shifted by - 0.1 w.r.t those of proton according 
to proton-to-helium ratio [5j. The Ferriere model [63J is used as the source distri¬ 
bution for the primary components, e.g. SNRs for SNe type II. To assure that the 
propagation parameters are correct, we need to compare the model prediction with 
the boron-to-carbon ratio and the proton flux [3]. As shown in Fig. [ 2 j the 
set of parameters used can reproduce the boron-to-carbon ratio and the proton flux 
well. According to this set of parameters, we can obtain the fluxes of the secondary 
positrons and electrons. A power-law spectrum with two breaks is introduced to pa¬ 
rameterize the injection spectrum of the primary electrons as a function of rigidity, 

f (p/Pbrl )" 71 (P<Pbrl) 

(p/Pbrl )" 72 (Pbrl < P < Pbr 2 ) ( 5 ) 

[ (Pbr 2 /Pbrl )' 72 ' (P/Pbr 2 )" 73 (P > Pbr 2 ) 

The parameters are adjusted according to the electron flux from AMS-02 [2], The 
agreement between the model and the data is shown in Section [5] These spectral 
indices are 71 = 1.95, 72 = 2.75 and 73 = 2.5 respectively. The breaks are pjj rl = 8.6 
GV and pjj r2 = 110 GV. Since the high energy breaks of primary particles, such as 
protons and helium, are found by PAMELA [ 6 ] and recently confirmed by AMS-02 
[3], it is reasonable to assume that there is also a high energy break in primary 
electron flux. More detailed discussion on the necessity of the high energy break 
p £ r2 can be found in [ 20 J and [ 68 ], where the high energy break hypothesizes are in 
favor compared to the no-break ones. Ref. [ 68 ] gave us an estimation by taking the 
primary electron flux as $ e - — $ e + and could roughly determine the break. 
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(a) 



Figure 2: (a) model prediction of B/C ratio compared with measurements from 
PAMELA [6], ATIC02 [7], CREAM-I [8] and TRACER06 [2j. (b) model prediction 
of proton flux compared with measurements from AMS-02 [3] . The solar modulation 
is taken as 500 MeV here. The red band in (a) shows the variation of the propagation 
parameters D 0 and 5 within 95% C.L. 


3 e ± from a single pulsar 

The pulsars are potential sources which could produce primary e 1 * 1 at high energy 
[33H361 Pf2| . Electrons can be accelerated by the strong magnetosphere of the pulsars, 
and this acceleration produces photons. When those photons annihilate with each 
other, they can produce e* pairs. Thus, the e ± energies are related to the pulsar 
magnetosphere. Assuming the pulsar magnetosphere as a magnetic dipole, this 
magnetic dipole radiation energy is proportional to the spin down luminosity. Due 
to this spin down (i.e. slowing of rotation), the rotational frequency of a pulsar 
Q = 2tt/P (with P being the period) is a function of time as follows [33, Ml 02j 

Do 

\/l + t/ r o 

where D 0 is the initial spin frequency of the pulsar and r 0 is a time scale which 
describes the spin-down luminosity decays. To cannot be directly obtained from 
pulsar timing observations, and it is assumed to be [351112] 

T 0 ~ 10 4 yr (7) 

The rotational energy of the pulsar is E(t) = (1/2)/D 2 (f). Here / is the moment of 
inertia, which is related to the mass and the radius of the pulsar and can be regarded 
as a time independent value. The magnetic dipole radiation energy is equal to the 



( 6 ) 





















energy loss rate, 


T 02 1 

l £ (‘)l = (8) 

The total energy loss of a pulsar is (35| SB, 02j 

£■.„.(*) = J‘df\E(t')\ = = |SWI« (l + )r) (9) 

The total energy injection of out of a pulsar should be proportional to the total 
energy loss 

E out (t) = rjE tot (t) = rj\E(t)\t fl + (10) 

where rj is the efficiency of the injected e ± energy converted from the magnetic dipole 
radiation energy. 

The pulsar characteristic age is defined as [58] 


(10) become 


p n 

— . — . — t T To 

2 P 2 Q 

(11) 

), we have T ~ t. In this condition, 

eqs. Q, (|9]) and 

\E(T)\ ~ 

(12) 

B„(T) ~ \E(T)\ — 

To 

(13) 

E„ C (T) ~ r,\E(T)\ — 

(14) 


To 


The propagation equation for the can be described as 


r) f f) 

-T = D(E)V 2 f + g^[b(E)f] + Q(x, t, E ), (15) 

where f(x,t,E ) is the number density per unit energy interval of e ± ; D(E) = 
{v / c)Dq{E /4 GeV) 5 is the diffusion coefficient with the velosity v of the parti¬ 
cle, the speed c of light, D 0 and S the same as the parameters used to calcu¬ 
late the background in Section [2j and b(E) = —dE/dt = b 0 E 2 with 6 0 — 1-4 x 
1CT 16 GeV _1 s _1 is the rate of energy loss due to inverse Compton scattering and 
synchrotron mmm- 

The source term Q(x,t,E) of a pulsar can be described by a burst-like source 
with a power-law energy spectrum and an exponential cutoff 


Q(x,t,E) 


QoE “exp 


E 

E cu t 



x 0 )S(t -1 0 ), 


(16) 
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where Q 0 is the normalization factor related to the total injected energy E out , a is 
the spectral index, and E cut is the cutoff energy. 

In Appendix [Aj we briefly review how to solve the equation (15) with the source 
([16]). The method is equivalent to many previous works (for example, Refs. [30; ^2]). 
Using the results, i.e. eqs. (56) and (60), in Appendix [A] we obtain the electron or 
positron flux observed at the earth as follows: 


^ c c Q 0 E 

$ e (r,t m ,E) = — f = - -s— 

47T 47T 7 T 2 r 3 , 


~ a ' E 


dif 


E n 


o-2 


exp 


E/E , 


cut 


d 2 


(1 — E / U max ) r 


dif 


(17) 

where d is the distance between the earth and the source, the diffusion distance r^a 
is given by 


r dif(hlif, E) — 2, 


D(E)t&i E ma 

\ (1 - 6) ~E 


E 


E„ 


i-5' 


(18) 


and the diffusion time f<jtf is the time a charged particle travels in the ISM before 
it reaches the earth. The electrons and positrons may be trapped in the pulsar 
wind nebula (PWN) for some time before they escape. The age of a pulsar is 
T = tescape + tdif, where t e scape is the time before the leptons escape from the PWN. 
In some case, Rscape and tdif can be of the same order of magnitude, and then 
the discussion will be complicated. In some other case, Escape could be negligible. 
For instance, when the SNR is evolving into the ” Sedov-Taylor” phase, the leptons 
in it are trapped (See Ref. [69] and references there in). In that case, the time 
^escape, during which the SNR reverse shock collides with the PWN forward shock, 
is typically a few 10 3 yr [@], which is small comparing to the ages of the pulsars we 
studied here, which are around 10 5 yr. In this work, we consider the latter case and 
neglect Escape for simplicity. We leave the case of large t esca , pe to a further specific 
study. Thus, we assume that tda — T. The maxium energy U max is defined as 


-Efnax — 1 /(bo T). 


(19) 


The positron fraction from AMS-02 implies a primary positron source with a cut¬ 
off energy 1 /E s = 1.84±0.58TeU _1 in their ’’minimal” model [I], which corresponds 
to E s G [490,790]. Due to the limitation of statistics of high energy e- and e+ 
measured by AMS-02, the upper bound 790 GeV is not a strict limit. Thus, we 
consider a primary e+ and e- source contribution with a cut-off energy U cut _ 0 g ~ 
(500 ~ 5000) GeV, which corresponds to a pulsar with an age T ~ (0.45 ~ 4.5) x 10 5 


yr according to (19). The term exp 


_ d 2 ' 

in ( 

17 

l r dif J 




in (17) tells us that a pulsar with d > r^f 


requires a larger normalization Q o, which hints a larger E out , a larger r/ in (10), or 
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both. r dif > d is required in our study, whose physical interpretation is that the 
distance a particle travels in the ISM should be larger than the distance between 
the earth and the source. Eq. (18) tells us that r d jf is as a function of diffusion time 
t d if and lepton energy E, as is shown by Fig. [3] where the color scale indicades r di f. 
For T ~ (0.45 ~ 4.5) x 10 5 yr and the lepton energy E = 1000 GeV, r di f is always 
greater than 0.5 kpc. Selecting pulsars with d < 0.5 kpc and T ~ (0.45 ~ 4.5) x 10 5 
yr, high engery leptons they produced can reach the earth. 



Figure 3: r dd as a function of i d jf and E, which is from eq. (18). The lepton energy E 
is the e + (or e~) energy detected at location away from the pulsar with the diffusion 
distance r dif . r dif increases with E. 


Thus, the pulsars with ages T ~ (0.45 ~ 4.5) x 10 5 yr and distance d < 0.5 kpc 
can explain the behavior of positron fraction of AMS-02 at high energy range. 


4 Single pulsar interpretation 

A few simple examples using a single pulsar are given to explain high energy positron 
fraction of AMS-02 [Tj. The background electrons and positrons are described in 
Section [2j The primary electron flux is scaled by a normalization factor A prirn e - 
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since it is not possible to constrain the electron flux contribution from SNRs. The 
age T and the distance d are taken from the ATNF catalogue and the positron 
fraction is fitted to obtain the free parameters in ( fl7| ) , the spectral index a, and the 
normalization Q 0 . Q 0 is fixed by the relation [27 ; , [36;, 41] E ont = /J™ dEEQ(E) ~ 
dEEQ(E), which approximately yields Q 0 ~ E out for a ~ 2. The cutoff energy 
E cut is set to be 5000 GeV, which is large enough, as it does not change the shape of 
pulsar contribution. Since we are interested in the positron excess at high energies, 
the fit is started from 10 GeV where the effect of solar modulation is negligible. 

Six nearby single pulsars, whose ages T ~ (0.45 ~ 4.5) x 10 5 yr and distance 
d < 0.5 kpc, are used to fit the positron fraction. Minuit package in ROOT is used 
to determine the parameters to minimize y 2 . The best results of the single pulsars 
are listed in Table [2j The results are also shown in Fig. [4} Using the parameters 


Pulsar name 

d(kpc) 

T(10 5 yr) 

O 

55 

o 

'--O 

a 

A 

Jr *- n pnm,e— 

X 2 /ndf 

Geminga 

0.25 

3.42 

50.5 

2.04 

0.50 

26.8/40 

J1741-2054 

0.25 

3.86 

50.6 

2.03 

0.50 

26.8/40 

Monogem 

0.28 

1.11 

50.1 

2.15 

0.50 

27.3/40 

J0942-5552 

0.30 

4.61 

50.6 

2.01 

0.49 

27.7/40 

J1001-5507 

0.30 

4.43 

50.1 

2.34 

0.47 

27.6/40 

J1825-0935 

0.30 

2.32 

50.5 

2.61 

0.44 

28.8/40 


Table 2: Parameters of six nearby single pulsars from the best fit results. The 
X 2 jndf from the fits of Geminga, J1741-2054, Monogem, J0942-5552 and J1001- 
5507 are smaller than 1, which show a good agreement between those single pulsar 
models and the experiment. 


of the best fit results, the positron fraction can be well reproduced by these single 
pulsar’s contributions. Table [2] tells us that the normalization A prirn e - are around 
0.5 and the spectral indices a of different pulsars are around 2. 

We can estimate the injection efficiency rj from the pulsar. Take Geminga as 
an example, the spin-down energy loss rate of Geminga \E(T)\ = 3.2 x 10 34 erg/s. 
The total radiation energy of the magnetic dipole can be derived from eq. (13) 
as E tot (T) ~ |F(T)|T 2 /r 0 = 1.2 x 10 49 erg. From the fit, we get the injection 
energy E out /2 = 10 5a5 GeV ~ 5.19 x 10 4 ' erg. From (14), we get r] ~ 8.7%. This 
efficiency is consistent with the previous studies by [33] and H2J. We can perform 
similar studies on the other five pulsars, whose results are listed in Table [3] A 
smaller i] means it is easier for this pulsar to produce the same amount of positrons 
and electrons. The efficiency required by J1001-5507 or J1825-0935 is too large to 
satisfy the physics condition for single pulsar interpretation. Geminga, J1741-2054, 
Monogem and J0942-5552 are the only candidates which survive from our selection 
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Figure 4: Single pulsar model can explain the positron fraction very well. According 
to the fitting result, the spectral indices are almost the same. 
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Pulsar name 

\E (10 33 erg) 

F tot (10 49 erg) 

£out(10 47 erg) 

v(%) 

Geminga 

32 

1.2 

5.19 

8.7 

J1741-2054 

9.5 

0.46 

5.83 

25 

Monogem 

38 

0.16 

1.90 

24 

J0942-5552 

3.1 

0.21 

7.02 

67 

J1001-5507 

0.68 

0.043 

1.89 

88 

J1825-0935 

4.6 

0.082 

5.09 

120 


Table 3: Electron injection efficiency r/ of the six nearby pulsars. For single pulsar 
interpretation of positron fraction, the results of Geminga, J1741-2054, Monogem 
and J0942-5552 are thought to be reasonable while the posibilities of J1001-5507 
and J1825-0935 as the high energy positron sources can be excluded. 


so far0 


5 Multiple pulsars interpretation 

The extra high energy positrons may come from serveral pulsars. We perform simi¬ 
lar study for multiple pulsars as we do for a single pulsar. Benefiting from the study 
in Section |4j, we can assume that the spectral indices a of all the pulsars are the 
same. Considering the physical models of the pulsars are similar, we make another 
assumption that the electron injection efficiencies rj are the same. These two as¬ 
sumptions help us reduce the number of free parameters. The discussion on rj from 
single pulsar in Section [4] tells us that Geminga, J1741-2054 and Monogem will give 
a much larger contribution to the high energy positron than J0942-5552. In other 
words, the rj of J0942-5552 in Table [3] is much larger than that of Geminga, which 
implies that the contribution from J0942-5552 in the multiple pulsars interpretation 
can be negligible compared with that from Geminga. 

We choose three from the four “surviving” pulsars in the multiple pulsars dis¬ 
cussion. The input parameters are the age T, the distance d and the energy loss 
rate E of each pulsar while the parameters we get from the fit is the normalization 
factor of primary electron A prim e -, the spectral index a and the electron injection 
efficiency rj. As shown in Fig. [5] (a), we obtain a good result from the multiple pulsar 
fit where y 2 /nd/' = 26.9/40. The parameters we get are A prim>e - = 0.50, a = 2.07 
and rj = 2.58%. The multiple pulsars interpretaion predicts a positron fraction with 

Considering that the uncertainty of log w {-Q^) from the fit is ±0.1, the E out for J1001-5507 
is 1.89/^36 x 10 47 erg . Thus, 77 = 88/^7 % for J1001-5507. There is no enough strong evidence 
that this 77 is smaller than 1. One should also note that 77 = 67/^ % for J0942-5552, which is 2cr 
smaller than 1 . 
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1 Electron Flux | 




Figure 5: Three pulsars fit to the positron fraction is presented here as an example 
of multiple pulsars fit, as show in (a). Using the same parameters, (b) shows the 
fitted parameters from positron fraction reproduce the electron flux when the solar 
modulation potential 550 MV is applied here. The error band in (a) shows the vari¬ 
ation of propagation parameters within 95% C.L., while that in (b) shows combined 
effect of propagation parameters and variation of solar modulation potential from 
400 MeV to 800 MeV. 


a decrease up to 600 GeV and after that a bump up to 2000 GeV, which is possible 
to be observed with more accumulating AMS-02 data. 

Using the parameters from the fit, we can reproduce the electron flux measured 
by AMS-02 [2] in Fig. [5] (b). It shows that our electron background estimation in 
Section 2 + pulsar contribution matches the experimental data especially at high 
energies. Fig. [5] also shows that the effect of the uncertainty due to the propagation 
model is small at high energy. The solar modulation potential is taken as 550 MV 
in the best fit result. The solar modulation potential is varied between 400 MV and 
800 MV to show that its effect on low energy is quite large. To reproduce the low 
energy electron flux more accurately, we need a monthly low energy electron fluxes, 
which may be published by AMS collaboration to model solar modulation. 

6 Discussion and conclusion 

In this work, we investigate the possibility that the rise of the positron fraction 
measured by AMS-02 can be explained by pulsars. The propagation parameters and 
the injection spectrums of nuclei and electrons are tuned according to the Boron-to- 
Carbon ratio and the proton flux. It will be better to tune those parameters with 
Boron-to-Carbon ratio and proton flux measured by AMS-02 since they are in the 
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same data taking period as the lepton fluxes. We find both the single pulsar model 
and the multiple pulsar model can explain the AMS-02 data very well. Six nearby 
pulsars are investigated as the single pulsar sources of the high energy positrons and 
finally four survive from all the conditions. The y 2 s of these single pulsars in this 
work are much smaller that those in [2T] . mainly because we set the cut-off energy 
equals to 5000 GeV while the authors of [21| set it to 1000 GeV. With three mostly 
contributing pulsars, the multiple pulsars model predicts a positron fraction with a 
decrease up to 600 GeV and a bump up to 2000 GeV. For the low energy, a simple 
solar modulation potential potential can not explain the measurement well. Thus, 
we need the monthly electron fluxes which can describe solar activity during the 
whole period. 

It is shown that the positron excess measured by AMS-02 can be explained by 
the pulsar scenario. Since the multiple pulsars can explain the experimental data 
well, it will be difficult to exclude pulsar scenario by isotropy. With accumulating 
AMS-02 data and future experiments, we can see the positron fraction behavior up 
to higher energy which will either confirm or reject the multiple pulsars scenario. 
If we consider other scenarios such as Dark Matter, we have to look into other 
productions, antiproton for instance, which have no contribution from pulsars. 


A Solving the diffusion energy-loss equation 


To fix the notation and for the pedagogical purpose, in this appendix we give a 
brief review on solving the diffusion energy-loss equation [30] j35] f38j [39l i42j- The 
diffusion energy-loss equation is given by 


' d_ 

Ft 


D(£)V 2 




S(x,t,E) = Q(x,t,E) 


( 20 ) 


where f(x , t, E ) is the particle number density per unit energy interval, D(E) > 0 is 
the diffusion coefficient, b(E) = —dE/dt > 0 is the energy loss rate, and Q(x,t,E) 
is the source term. 


A.l Green function for the diffusion energy-loss equation 


The Green function G(x,t, E; x 0 ,t 0 , E 0 ) of (20) is defined as 


' d_ 

Ft 


D(E)V 2 




G(x,t,E-,x 0 ,t 0 ,E 0 ) 


5 3 (x-x 0 )5(t-t 0 )5(E-E 0 )(21) 
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The solution of (21) has been given in Ref. [BB]. To fix the notation, let us briefly 
review the derivation. Define 


cj)(x, t, E] x 0 , t 0 , E 0 ) = b(E)G(x, t, E; x 0 , t 0 , E 0 ) 


( 22 ) 


Substituting G = (j>/b into (21) gives 


d_ 

dt 


b ^WE 


D(E)V 2 


4>{x, t, E- x 0 , t 0 , E 0 ) = b(E)5 3 (x 


x 0 )5(t - t 0 )S(E 


Let us make the variable transformation (t, E) —y (t/, A) as follows, 


-E 0 ) 

(23) 


t' = t — r(E, E 0 ) , with t{E, E 0 ) = 


cEq 


dE 

W) 


X(E,E 0 ) = 


cEq 


dE 


D(E) 

b(E) 


The Jacobian matrix of this transformation is easily obtained as 
d(t', A) 


whose inverse is 


dit, E) 




EL 

dt' 

dt 

dX 

dE 

dX 

dt 

dE 

dt 

dt 

dt' 

dE 

dX 

dE 

dt' 

dX 


l 


b (E) 

D(E) 

b(E) 


Thus, 


d(t\ A) - Vf t 


d dt d dE d 


0 


D(E) 
b(E ) 
D(E), 


d A dXdt dX dE D(E) [dt 


d 


~ b{E) m 


Substituting the relation d t — b{E)dE = D(E)d\ into (23) implies 

-<5 3 (f - f 0 )6(t - t 0 )5(E - So) 


d 2 
-V 2 

g\ 


KE) [3 , 


D(E) 
b(Eo ) s3 


D(E„ 


6 J (x - x 0 )8(t - t 0 )S(E - S 0 ) 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 


It follows from eq. (26) that 


det 


~d(t', xy 

d{t, E)_ 


D(E) 

m 


(30) 
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which implies 


5(t-t 0 )5{E-E 0 ) = 


D(E 0 


b(E 0 ) 
D(E 0 ' 
b(E 0 ) 


bit' -1 0 )5(\) 
S(t' - t 0 )5(X) 


(31) 


where in the second equality we have used the properties: D(E 0 ) > 0 and b(E 0 ) > 0. 


Substituting (31) into (29), we obtain 


f - V 2 i <j> = S(f - to) S 3 (x - x„)S(\) 


(32) 


As is well known, the Green function, which can be written as G[x — x o, A) in 
the spherically asymmetric approximation, of the diffusion equation satisfies 


- V 2 ) G(x - x 0 , A) = <5 3 (f - x 0 )S(A) 


with A > 0. The solution of (33) is 


G(x — x 0 , A) = 


(4ttA) 3 / 2 


exp 


(x - Xo f 
4A 


(33) 


(34) 


Comparing (32) with (33), we can read off the solution of <fi as follows 

{x - x 0 ) 21 


= S(t'-t 0 ) 


(47tA) 3 / 2 


exp 


= 8[t - t 0 - t(E, E 0 )} 


(47tA) 3 / 2 


4A 

exp 


(x-xp f 
4A 


(35) 


where in the second equality we have used (|24|). Thus, we finally get the solution of 

(x - x 0 ) 2 


(21) as follows 

G(x,t,E-,xo,t 0 ,Eo) = 


4> _ 5[t - t 0 - r(E,E 0 )} 

b(E) ~ KB)[4xA(B,£: 0 )]3/ jeXP L 4A(B,B 0 )J 


(36) 


where the functions t(E,E 0 ) and A{E,E 0 ) are defined in eqs. (24) and (25), respec¬ 
tively. 

A.2 Solution for a burst-like source 


Once we know the Green function, eq. (|36|), we can write down the solution of 

(37) 


eq. (20) for a generic source as follows 

f(x,t,E)— [ d 3 x'dt'dE'G(x, t, E-, x', t', E')Q(x', t', E') 
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In particular, for a burst-like source, the source function is proportional to S 3 (x — 
xo)S(t — t 0 ), that is, 


Q(x, t, E) = Q(E)5 3 (x - x 0 )5(t - t 0 ) (38) 

where Q(E) is an arbitrary function of E, x 0 is the position of the source, and t 0 is 
the instantaneous time when the source bursts. Substituting eqs. (36) and (38) into 
eq. (37), we obtain 

E) 


b(E)[47rX(E,E')}y 2 


f ^,5[t-t Q -r(E 1 E')} 

= I dE . , r . — — exp 


exp 


(x — x 


A 2 


4A (E,E') 


(x-x 0 f 


&(£)[47tA(£,£')] 3/2 L 4A (E,E')_ 

Denote the solution of the equation 

0 = t — t 0 — t(E, E') = t — t 0 — 


Q(E') 


Q(E')5 3 (x' - x 0 )6(t' - t 0 ) 


(39) 


fE' 


dE 

W) 


is E 1 — Eq, then we have 

S[t - t 0 - r(E, E')} = b(E 0 )5(E' - E 0 ) 

where we have used 
dr 


dE' 


5[t-t 0 -r(E J E')\ = 5(E' -E 0 ) 


(40) 


(41) 


(42) 


E'=E 0 


Substituting eq. (41) into eq. (39), we have 

^ f b(E 0 )8(E'- E 0 ] 
f (x, t, E) = / dE - , exp 


b{E)[47T\(E,E')]y 2 4A(£,£') 


(x- x 0 y 


Q(E') 


Q(E 0 ) b(E 0 ) 

[AeX(E,E 0 )}V 2 b(E) 


exp 


(X - Xp) 

’AX(E,Eq) 


2 1 


(43) 


where the initial energy Eq is defined as the solution of eq. (J40|). In other words, if 

(44) 


we know t — t 0 and E, we can find E 0 by solving the equation 

r Ro dE 


t - t 0 = t(E, Eq) = 
Define the diffusion distance r^if as 


m 


"Eq 


r dii (E, Eq) = V4A(£,£ 0 ) = 2\ dE 


D(E) 

b(E) 


(45) 
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we can rewrite eq. (43) as 


f(x, t, E) = 


QCgo) b(E 0 ) 

7 rir^ f b(E) 


exp 


dif 


(46) 


with r 2 = (x — x o y 


A.3 Solution for a burst-like source with power-law spec¬ 
trum 


Consider the case when the function Q{E ) in eq. (38) is a power-law function, that 
is, 


Q(E)=Q 0 E~ a ex p 


E 


E, 


cut 


(47) 


where Q 0 is a normalization constant, a is the spectral index, and E cut is the cutoff 
energy. The diffusion coefficient D(E) and the energy loss rate b(E) are assumed to 
take the form 

(48) 

(49) 


D(E) = PD„ (T 
b(E) = b 0 E 2 


where j3 = v/c is the ratio of velocity to speed of light, the constant D 0 and the 
index 5 can be figured out by the background fitting in Sec. [2j the constant b 0 is 
given in Sec. [3j and E 0 should be determined by eq. (44). 

Let us calculate the initial energy Eq first. Denote the diffusion time tdif = t —to- 


It follows from eqs. (44) and (49) that 

r Eo dE 


tdif — 




E Er 


which implies 


i “ w 0 - iotm 


(50) 


(51) 


from which, we see ^ > b 0 tda, that is, E < Denote E max = l/(6 0 tdif), then 

= (52) 
u U/o Dnax 

which gives 


E n = 


EE n 


f'max E 


(53) 
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Thus, we obtain the following pieces in eq. (46): 


Q(Eq) = QqE 0 a exp 


= Q 0 E~ a 1- 


Eq \ _ „ ( EE max 

— Vo 


E, 

E 


cut 


E, 


E n 


exp 


- E 

max 

E / E cut 


exp 


EE n 


E cut (E max E) 


(1 — E/E max ) 


( 54 ) 


2 / E 


b(E 0 ) E 0 
b(E ) E 2 


E u 


-2 


(55) 


Substituting eqs. (54) and (|55j) into eq. (46), we obtain 

Q 0 E~ a f t e v 2 r E / E 


f(x, t, E) = 


7 T2 r 


dif 


E n 


exp 


_ (1 E / E m ax ) 


(56) 


which is consistent with previous works (for example, eq. (6) of Ref. [42]). 

Now let us figure out the diffusion distance r dif . To this end, we substitute 


eqs. (48) and (49) into eq. (25), and get 


A (E,E 0 ) = r ^ r d EE‘-> 

y ] Je b 0 E 2 Efoo J E 

PDo / 77,5-1 77.5-1A 


E S M5 - 1 ) 


(Et 1 ~ E 6 - 1 ) 


(57) 


which can also be written as 


\(e,e 0 ) = pd 0 [E) 


7-dif 


Eq J b 0 tda (S — 1) E 




(58) 


Comparing the above equation with eq. ( |48| ), E max = 1 /(&o^dif) and E 0 /E = (1 — 
E/E max ) _1 , we obtain 


KE,E 0 ) = 


D(E)tE d 


(1-5) E 
Thus, the diffusion distance r d if is given by 


1 - 1 - 


E 


E„ 


1—5 


(59) 


T, dif = 2 a/A (E, Eq) — 


D(E)tdii E max 

1 

1 

i_ 

(1-5) E 

\ ^max / 


(60) 


which is consistent with previous works (for example, eq. (10) of Ref. [33] and eq. (7) 
of Ref. [42]). 
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